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ABSTRACT 



A mathematical model is developed to analyze the seismic 
response of a drydocked vessel in three degrees of freedom; 
vessel rotation about its keel and vessel horizontal and 
vertical translations relative to the drydock cradle. Data 
from eleven actual vessel-drydock systems and the time 
acceleration history of an earthquake are implemented to 
predict vessel three degree of freedom response during an 
earthquake. Vessel seismic response and the resultant drydock 
forces are compared to vessel-drydock system failure criteria 
to determine stability of the system during earthquakes. The 
three degree of freedom vessel response model is compared to a 
one and two degree of freedom vessel response models and a 
model in which seismic loading is simulated by a single static 
force. The three degree of freedom vessel motion is shown to 
be the most accurate method for analyzing vessel-drydock 
system failure criteria. 
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1.0 INTRODUCTION 



1.1. Background 

In the design of a drydock for military vessels, an 
important factor is the environmental-generated forces that 
are found in the drydock cradle which supports the ship. 
Earthquakes and the wind are considered the two most 
significant environmental loading factors. This thesis 
explores the seismic loading aspect of the problem. 

Regardless of its location, any U.S. Navy ship that is 
nuclear powered and in drydock is considered to be in a high 
seismic risk area. Currently the Navy defines a high seismic 
risk area as one in which earthquake forces be approximated by 
a steady horizontal force of 0.2 g times the vessel mass 
acting at the center of gravity. The drydock can be analyzed 
to see if it can withstand this quasi-static loading. 

1.2 Previous Work 

In 1981 B.V. Viscomi studied the seismic response of a 
drydocked submarine using a single degree of freedom model 
[14]. This analysis assumed that the vessel was allowed to 
rotate about the keel. Viscomi used the time history 
acceleration record of the North-South component of the 1946 
El Centro, CA Earthquake, California Institute of Technology 
processing scheme. The vessel was analyzed for system failure 
(i.e. the vessel lifting off of a row of side blocks). 

A thesis by C.F. Barker in 1985 [1] used a two degree of 
freedom model to study seismic response. Employing the 1946 
El Centro Earthquake, Massachusetts Institute of Technology 
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standard processing scheme, Barker analyzed rotation about the 
keel plus horizontal relative translation between the vessel 
and the drydock cradle perpendicular to the vessel longitudal 
axis. The two degree of freedom model permitted the analysis 
of drydock cradle failure: block sliding and tipping, which 

were not possible to analyze in Viscomi's one degree of 
freedom model, along with block liftoff. 

1.3 Contributions of This Thesis 

This thesis develops a three degree of freedom model of a 
vessel in drydock. Motion will be permitted in three 
directions: rotation about the keel and relative 

translations, horizontal and vertical, between the vessel and 
the drydock cradle. The three degree of freedom model will 
make possible the analysis of drydock cradle failure: 
additional block crushing forces due to vertical motion, which 
was not conceivable in the previous one and two degree of 
freedom models developed by Viscomi and Barker respectively, 
along with block sliding, tipping and liftoff. 

The seismic input for the three degree of freedom model 
is the time acceleration history of the North-South component 
of the 1946 El Centro, CA Earthquake, Massachusetts Institute 
of Technology standard processing scheme. The MIT version of 
the El Centro Earthquake acceleration record is described in 
reference [9] , and is displayed graphically in Figure 1.1. 

1.4 Outline of This Thesis 

Section 2 of this thesis describes the vessel-drydock 
system for the reader and introduces terms which are used 
throughout the thesis. Section 3 describes the various 
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failure mechanisms found in the vessel-drydock system which 
can occur. 

Section 4 examines the modelling of the vessel-drydock 
system under seismic loading. Section 4.1 describes the 
quasi-static loading method, and Section 4.2 develops the one 
degree of freedom equation of vessel motion model. Section 
4 . 3 develops the two degree of freedom equations of vessel 
motion model. Finally, Section 4.4 develops the three degree 
of freedom equations of motion model. 

Section 5 discusses the system parameters that are 
required for the mathematical model found in Section 4. 
Included in this Section is the modelling of block stiffness. 
Eleven typical vessel-drydock system configurations are 
studied to be implemented into Section 6. 

Section 6 examines methods of evaluating the one, two and 
three degree of freedom models found in Section 4. Section 
6.1 evaluates the linear equations of motion models using 
modal analysis method. Section 6.2 evaluates both linear and 
non-linear equations of motion models using a fourth-order 
Runge-Kutta numerical analysis scheme. Section 6.3 explains 
the response spectrum method of determining the maximum value 
of vessel seismic response. Section 6.4 determines system 
response using quasi-static force method. 

Section 7 describes the development and testing of the 
one, two and three degree of freedom vessel motion computer 
programs, and contains the results for a seismic analysis of 
several actual vessel-drydock systems. 

Section 8 summarizes the response predictions for several 
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actual vessel-drydock systems obtained with quasi-static, one 
degree of freedom, two degree of freedom and three degree of 
freedom analysis methods, and conclusions are drawn. 

Appendix 1 is an example listing of the fourth-order 
Runge-Kutta computer program. Appendix 2 is an example 
listing of the modal analysis computer program. Finally, 
Appendix 3 describes the modal analysis of the two and three 
degree of freedom vessel-drydock systems and predicts maximum 
system response using the response spectrum method of analysis 
and modal participation factors. 
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FIGURE 1.1 

Acceleration Time History, El Centro Earthquake, 1946 
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2.0 THE VESSEL-DRYDOCK SYSTEM 



2.1 General Description 

In order to perform inspection, maintentance and repair 
on the outer hull of a ship, it is necessary to place the 
vessel in a drydock. A drydock is a concrete encasement built 
into the earth. At one end of the drydock is a flood gate to 
allow the entry and departure of the vessel. A cradle, formed 
by drydock blocks, is built upon the drydock floor in order to 
support the ship once the water is pumped out of the 
encasement. The drydock must be able to hold the weight of 
the vessel and cradle. 

Due to the various hull configurations associated with 
the different classes of ships, the cradle component of the 
drydock-cradle-vessel system is varied to suit a particular 
docking situation. The cradle accomplishes the transfer of 
vessel weight to the drydock floor, provides stable support to 
the vessel and allows access to the vessel hull. This cradle 
is usually constructed of timber, concrete, or a timber- 
concrete composite. 

2 . 2 System Component 

The primary components of the drydock-cradle-vessel 
system are shown in Figure 2.1. This section defines terms, 
related to docking, used in this thesis to describe the system 
components . 
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1. Blocks - The units, consisting of timber, concrete, steel 
and other materials, which together make up the cradle. 

2 . Cradle - A framework of blocks which supports the vessel 
when the drydock is dewatered. 

3. KG - The distance between the vessel baseline (ship's 
keel) to its vertical center of gravity. 

4. Keel Blocks - The center blocks, directly beneath the 
vessel's keel. 

5. Pier - A column built of blocks that extends from the 
ship hull to dock floor. 

6. Side Blocks - The blocks located to the right and left of 
the keel blocks. 

7. Ton - A long ton, 2240 lbs. 

8. g- Acceleration of gravity, 32.2 ft/sec 2 . 
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3.0 VESSEL-DRYDOCK SYSTEM FAILURE MECHANISMS 



There are four failure inodes that the vessel-drydock 
system can exhibit: block crushing, block sliding, block 

overturning and block liftoff (i.e. the liftoff of vessel from 
a row of blocks) . 

3 . 1 Block Crushing 6 

A block can support compressive stresses linearly until 
the blocks' proportional limit is exceeded. The proportional 
limit is a material property of the block and is the maximum 
compressive stress at which stress is still linearly 
proportional to strain. If the compressive proportional limit 
is exceeded, the block is considered to have failed. 

3.2 Block Sliding 6 

Due to friction at the interfaces between the blocks and 
the vessel and between blocks with other blocks, drydock 
blocks will have the tendency to resist sliding when subjected 
to vertical and horizontal loads. This is true when in the 
absence of mechanical fasteners. 

First consider the blocks that form the keel pier. 

Figure 3 . 1 shows that the resistance of the keel pier to 
sliding is 

H = y 2.V (block-block interface) 
and H = y 2 V (vessel-block interface) 
where 

H = horizontal resisting force 

= coefficient of friction, block-block interface 
y 2 = coefficient of friction, block-vessel interface 

V = vertical load on the blocks. 
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Failure of the keel pier due to sliding will occur when 
the ratio of horizontal force to the vertical force is greater 
than the lowest coefficient of friction found in the block. In 
other words, sliding occurs when 
~ >(y 1 or y 2 ) . 

As with the keel piers, the resistance of the side piers 
to sliding will depend on the interface which offers the least 
resistance, the block-vessel interface or the block-block 
interface. However, this case is more complicated due to the 
complex geometry of the side pier as shown in Figure 3.2. 

The block-vessel interface of the side pier is examined 
and shown in Figure 3.3. An arbitrary force F at some angle 
from vertical is applied to the face for the side block cap. 

At this interface, the horizontal and vertical block force 
reactions are: 

Normal Force = H sin4> + V cos4> 

Tangential Force = H cos<{> - V sin4> 
where H = Horizontal Reaction 

V = Vertical Reaction 
$ = Block Inclination Angle. 

Using the formula for horizontal sliding resistance, H = y 
V, previously introduced in the keel pier case and the 
above equations yields 

H cos<})'- V sintp = V 2 (H sin<}) + V coscj>) . 

Rearranging the above equation gives 

H _ P2 cosd> + sin<t> 

V costf - U 2 sin<t> 

or 
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where 



y 2 cos<J) + sin<l> 

^2 ~ cos^ - P2 s; '- n ^ 

For the block-block interface of the side pier, sliding 
resistance is found by taking the vertical and horizontal 
components of the applied force, F, and using them in the 
formula H = y V. Like in the keel pier, failure of the side 
pier due to sliding will occur when the following is satisfied 
^ or y 2 ' ) • 

There is a critical angle a where no slippage occurs. 
Referring again to Figure 3.3, this angle can be calculated by 
comparing the maximum tangential force F^ max the interface 
can have without slipping and the corresponding normal force 
F n . The following relationship holds 
Ft max = ^F n 



Finally, 

a = arc tan y. 

Since the applied force F acts in a straight line through all 
blocks to the ground, the angle oC for each interface must be 
calculated. If force F is applied outside of these angles, 
slippage will occur. 

3.3 Block Overturninq 11 

The third failure mechanism is the overturning of a block 
due to an applied force. The line of action of this force 
must fall within the middle one-third of the base of the 
block, as shown in Figure 3.4, or the block will tip over. 

In order to show the region that the block will remain 
upright, the inclined force F must be broken into a transverse 
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component and a vertical component F v . Superposing the 
bending and axial stresses caused by F^ and F v gives the 
resultant stress at any point in the block as 




where L = block height 

A = base area, here bh 
I 2 = base moment of intertia, here ^ 3 /12. 

The minimum stress point on the bottom of the block 
occurs at point M (Figure 3.4), a distance -h/2 from the z 
axis. Hence, 

/_ x _ v 

°x m bh h bh 

When ° XItl at M is negative, compressive stress is present and 
the block is stable. When at M is positive, tensile 
stress is present and the block will overturn since no 
fasteners hold block to the ground to develop tensile stress. 
The limit of block stability is when 

(°x)m = 0 • 

Thus, the following equation holds 

F = ^ F . 
h 6L v 

This condition exists when force F is applied at an angle such 
that the line of force lies within the one-third of the base 
of the block. If the line of force is outside this region, 
the block will tip. 

3.4 Vessel Liftoff 

The fourth and final failure mode occurs when the vessel 
breaks contact with either the side or keel piers. This 
failure occurs when the dynamic deflection of a row of blocks 
is equal to, or exceeds the average static deflection of the 
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blocks. The static deflection of blocks is given by 



W 



6 s 2K + K 
sv 



kv 



where W = submarine weight 

K sv = side pier vertical stiffness 
K kv = keel pier vertical stiffness * 



The dynamic deflection can be due to vertical displacement y, 
rotational displacement 6, or a combination of both. 
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FIGURE 3.3 

Side Pier Forces, Block-Vessel Interface. 



- 25 - 




FIGURE 3.4 

Block Tipping Forces 
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4.0 MATHEMATICAL MODELS OF VESSEL-DRYDOCK SYSTEM 



This portion of the thesis discusses the modelling of the 
drydocked submarine subjected to a seismic load. Models will 
include a quasi-static method, one, two and three degree of 
freedom mathematical models. The models developed in this 
section will be analyzed in Section 6 and 7 to predict the 
vessel seismic response. 

4.1 The Quasi-static Method of Modelling Seismic Response 

4.1.1 Approach 5 < 6 

Current U.S. Navy design method for submarine's seismic 
response will be examined in this section. The submarine is 
described as a rigid cylindrical body with its weight evenly 
dispersed longitudinally. 

The quasi-static force method replaces the earthquake 
motions by a force corresponding to the vessel mass time 
0.2 g. This force is horizontally applied to the submarine's 
center of gravity in the tranverse direction. The drydock 
blocking system is determined by the evaluation of the loads 
generated by the quasi-static force. 

4.1.2 Force Equations 

The application of the static force to the center of 
gravity of the submarine in drydock is represented in Figure 
4.1. The seismic overturning moment, M s , is defined as 

M s = (A/g) (a) (KG) (2240) ft-lb 

where 

A = vessel displacement in long tons 
g = acceleration of gravity 

a = vessel's center of gravity acceleration. 
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Current U.S. Navy design practice states that the acceleration 
of the submarine's center of gravity, a, should be set equal 
to 0.2 g. Hence 

M s = 448(A) (KG) ft-lb 

After determining the seismic overturning moment, the 

number of sideblocks required on one side of the vessel is 

M 

s 



where N = # of side blocks required in the row 

A = contact area of a block in square inches 

Sp = proportional limit of the block cap in lb/in 2 

L = distance between centerline of keel and side 
blocks, in feet, as shown in Figure 4.1. 

4.1.3 Vessel Response in the Quasi-static Load 
Method 

The application of a force applied at a vessel's center 
of gravity will generate reactionary side block forces in the 
vertical direction to oppose the seismic moment, as shown in 
Figure 4.2. Summation of moments about the keel yields the 
equation 

ZM = M s - L(F sr + F mr ) + L(F sl - F ml ) = 0 (4.1.1) 



where 

F s = static pier forces (left and right) 

F m = pier forces due to the applied moment (left and right). 
Due to symmetry in the vessel-drydock system, the 

applied moment resistive forces in the side blocks are equal 
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in magnitude but opposite in direction. Thus, Equation 4.1.1 
simplifies to the form 



M s = 2L(F m ) (4.1.2) 

where F m = F m ^ = F mr . 

The force F m is equal to the side pier compression 
displacement, 6m, times its vertical stiffness, K vs . This 
compression displacement can be expressed as 

6 m =Lsin(G). (4.1.3) 

where 0 = angle or vessel rotation about its keel. 

Assuming the rotational angle is small, Equations 4.2.1 and 
4.1.3 are combined to yield 



M 



e = 



2 

2L K 



sv (4.1.4) 

4.2 One Degree of Freedom: Rotational Response 

4.2.1 Approach 

A one degree of freedom mathematical model of a submarine 
in drydock is developed in this section. The submarine is 
identified as a rigid cylindrical body with an even 
longitudinal weight distribution. The system damping and 
stiffness are provided by the drydock blocking arrangement. 
System excitation is due to seismic ground accelerations. 

Other assumptions are 



The keel and side piers remain vertical during the 
ground motions, 

no slippage between the block cradle and drydock, and 

the keel and side piers bases accelerate at the same 
rate as the drydock (ground) . 

The idealized model of the vessel-drydock system is shown 
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in Figure 4.3. In this system, vessel motion is restricted to 
rotation about the keel. The side piers are modelled as 
vertical due to typically small angles or inclination of the 
tops of the side piers. Only vertical displacements in the 
side piers are considered. 

The one degree of freedom mathematical model predicts and 
analyzes vessel rotational response about its keel caused by 
seismic ground accelerations normal to the submarine's 
longitudinal axis. The resultant angle or rotation, expressed 
in radians, and side pier forces, generated by the rotation, 
are compared to applicable pier failure criteria found in 
Section 3, to see if the blocking arrangment will remain 
intacted. 

4.2.2 One Degree of Freedom Equation of Motion Model 
For the one degree of freedom case, the keel is the 
origin with clockwise rotation, 6 , positive. The simplified 
vessel-block system showing coordinate 6 at rest and when 
excited are shown in Figures 4.3 and 4.4 respectively. The 
moment equilibrium equation about the keel (origin) is: 

EM k = I k 0 + MKGXg - MKGy g sin6 

= B/2 (F lsv - F ldv ) - B/2 (F 2sv + F 2 dv) + WKG sin 0 (4.2.1) 

where I k = Mass moment of inertia of vessel about the keel 
M = Vessel Mass 

F lsv' F 2sv = static side block forces in the vertical direction 
(Figure 4.4) 

F ldv' F 2dv = dynamic side block forces in the vertical 
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direction (Figure 4.4) 

W = vessel weight (Mg) 

and Xg,yg = seismic ground accelerations in their respective 
directions 

This equation can be simplified by the following: 

1. ) The left and right side block rows have the same 

vertical spring constant, K sv , have the same static 
deflection and equal but opposite dynamic deflection due 
to symmetry, i.e. 

^isv = F 2sv = F sv (4.2.2) 

and I F ldvl = l F 2dvl = F dv* (4.2.3) 

2. ) The dynamic side block force, F<} v , as seen from Figure 

4.5 equals to a force due to the modelled spring 
displacement plus a dissipative force due to system 
damping. In other words, 

F dv( s P rin 9) = (B/2)K SV sin6 
and F^ v (dissipative) = C(B/2) 6 
or when combined, 

Fdv = B / 2 ( K SV sin6 + c6 ) (4.2.4) 

where K sv = vertical spring constant of the side blocks 

C = vertical damping coefficient of the side blocks 

Substituting Equations 4.2.2, 4.2.3, and 4.2.4 into 

Equation 4.2.1 and rearranging components yields 
•• • 

1*6 + ( B 2 / 2 ) C 9 + [(B 2 /2)K sv - W KG] sin6 

= -MKG X g + MKG y g sin 0 (4.2.5) 

An additional simplification is that since 0 is very small, 
sin 0 ~ 0. 

The final simplification is that damping is expressed as a 
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fraction of critical damping. To do this, allow 

(B 2 /2 ) C =5 C cr = 2? I k oh 

where £ = viscous damping factor(B 2 C/2C cr ) 

C cr = critical damping coefficient = 2I k w n 

w n = system undamped natural frequency. 

Now, Equation 4.2.5 simplified to 

I]c S + 2 5 I K “ n 9 + [ (B 2 /2 ) K sv - WKG]0 

= -MKG X g + MKG y g 0 (4.2.6) 

The non-linear term, MKG y g 0, found in the equation of 

motion will be removed at this time to linearize Equation 

4.2.6 since 6 is assumed to be small. The effects of the 

non-linear term will be evaluated in Section 6 and 7. 

The final form of the linearized one degree of freedom 

equation of motion is 
• • • 

m l® + c l 9 + k i© = -m 2 Xg (4.2.7) 

where m;L = I k 

C 1 = 2 ^ ^k^n 

k! = (B 2 /2 )K sv - W KG 

rn 2 = M KG. 

In matrix form, 




To simplify the explanation, the matrices will be redefined as 



[A]{y) + [B] (y) = (E(t) } 



(4.2.8) 



where 

[A] 
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• 

0 




~ o 


{y} - 


6 


{E (t) } = 


-m n x 








2 g 



Equations 4.2.6 and 4.2.8 will be solved in Sections 6 and 7. 

4.3 Two Degree of Freedom: Rotational and Horizontal 

Translation Response 

4.3.1 Approach 

This section develops a two degree of freedom 
mathematical model of a drydocked submarine. As in the case 
of the one degree of freedom (Section 4.2.1), the model 
assumptions still hold, i.e. a rigid cylindrical body, keel 
and side piers remain vertical, no slippage between the block 
cradle and drydock, and the bases of the keel and side piers 
accelerate at the same rate as the drydock (ground) . 

The idealized model of the vessel-drydock system is shown 
in Figure 4.6. In this system, vessel motion is restricted 
to rotation about the keel, and in horizontal transverse 
translation relative to the keel and side pier supports. The 
approach for determining vessel-drydock failue is identical to 
the one degree of freedom case except for the addition of 
sliding and tipping failure modes. 

The two degree of freedom mathematical model predicts and 
analyzes submarine rotational and transitional response caused 
by seismic ground accelerations normal to the vessel's 
longitudinal axis. The rotational response is in radians, the 
translational response in inches, and both responses are 
expressed as a function of time. The keel and side pier 
forces are compared to applicable pier failure criteria found 



- 33 - 



in Section 3. 



4.3.2 Two Degree of Freedom Equation of Motion Model 
As in the one degree of freedom case, the keel is the 
origin for the two degree of freedom system with clockwise 
rotation, 6 , is positive. The linear relative translation 
coordinate, x, is defined with respect to vessel and drydock 
translation relative to the ground. Defining u as the 
position of the submarine keel relative to the ground, then 
the following relationships hold: 

x = u - x g 

U = X + Xg 

and u=x+Xg. (4.3.1) 

The simplified vessel-block system showing coordinates (6 and 
x) at rest and when excited are shown in Figures 4 . 6 and 4 . 7 
respectively . 

A balance of forces in the x direction, as shown in 
Figure 4.7, yields the first equation of motion 

IF X = Mu + MKG 0 = (F lsh - F ldh ) - (F 2sh + F 2dh ) 

- (F 3 sh + F 3 dh) (4.3.2) 

where F-|«-h/ F 2 eh/ F^-h = static block forces in the horizontal 

direction at their respective 
position 

Firih/ Forthi Fthh = dynamic block forces in the 

horizontal direction at their 
respective position. 

The above equation can be simplified by the following: 

1) The left and right side block rows have the same 
horizontal spring constant, K S h, have equal but opposite 
static horizontal deflection and the same dynamic horizontal 



- 34 - 



deflection due to symmetry, i.e., 

F lsh + F 2sh = 0 

F ldh = F 2dh • (4.3.3) 

Also, the keel block rows do not experience any net static 
horizontal deflection since the inclination angle of the block 
cap is zero. Hence 

F 3sh = 0* (4.3.4) 

2) The dynamic block forces, F^ as seen from Figure 4.7 
equals to a force due to the modelled springs plus a 
dissipative force due to system damping. In other words 

F dh (spring) = F ldh (spring) + F 2dh (spring) + F 3dh (spring) 
= K sh x + K sh x + F kh x 

F dh (spring) = ( 2K sh + K kh) x 
Likewise 

F dh (dissipative) = c x x 
so 

F dh = ( 2K hs + K hk) x + ^x x (4.3.5) 

where Kg^ = side pier stiffness 

K kh = ^ ee l pier stiffness 

C x = system horizontal damping coefficient. 
Substituting Equations 4.3.1, 4.3.3, 4.3.4, and 4.3.5 
into Equation 4.3.2 and rearranging components yields 

mu x + m 12 ® + Ci x + kn x = -mu x g (4.3.6) 

where m^ = M 

mi 2 = M KG 

C 1 = c x = system horizontal damping coefficient 
kn = (2K s i 1 + Kjth) = system horizontal stiffness 
Summing the moments about the origin in a similar fashion 
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as the previous degree of freedom case, as shown in Figure 
4.7, yields the second equation of motion, 

IM K = I k 0 + MKG u - MKG y g sin0 

= B/2(F lsv - F ldv ) - B/2 (F 2sv + F 2dv ) + WKG sinQ (4.3.7) 

The second equation can be reduced as it was in the one degree 
of freedom case to 

I k e + MKG X + C e ® + (B 2 /2 K sv - WKG) 6 

= MKG X g + MKG y g 6. (4.3.8) 

Once again, the non-linear term, MKG y g 0, found in 
Equation 4.3.8 will be removed leaving a linearized equation 
since 0 is assumed to be small. The effects of the nonlinear 
term will be evaluated in Sections 6 and 7. 

Now the system of equations for the two degree of freedom 
vessel-drydock system is as follows, 

mu x + m 12 9 + Cn x + kn x = “^ii * g (4.3.9a) 

and 

* * • • • • • 
m 22 6 + m 21 x + c 12 6 + k 22 0 = ~™21 x g (4.3.9b) 

where = M 

m 2 i = irii2 = m KG 
m 22 = I k 

c± = system horizontal damping coefficient 
c 2 = system rotational damping coefficient 

k ll = 2K sh + K kh 
k 22 = B 2 /2 K sv - WKG. 

These two equations of motion are coupled in the mass times 
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acceleration term, which is known as inertial coupling. Also 
the equations are coupled in the damping term, because the 
system's two natural modes have translation and rotation 
involved in them. This model uses 5 percent of the critical 
damping coefficient as system damping. To evaluate damping, a 
modal analysis for the system is performed. This analysis is 
performed in Appendix 3. Substituting damping coefficients 
into Equations 4.3.9a and 4.3.9b, the equations of motion can 
be expressed in matrix notation as Equation 4.3.10: 



[M] (yJL) + [C] (y_L) + [K] (yl) = (E'(t)} (4.3.10) 

\ 



where 



(y'} = Response Vector 




[M] = Mass Matrix = 



m ll ^12 
L m 21 m 22. 



[C] 

[K] 



- Damping Matrix = 



= Stiffness Matrix 



’ C ll C 12 _ 
- C 21 C 22 . 




{ E ' ( t ) } = Seismic Forcing Vector 





To simplify the explanation, the matrices will be redefined as 
[A] {y} + [!] {y} = { E ( t) } (4.3.11) 

where 



[A] 



o [M] " 
_[M] [C] _ 
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- [M] 



O 



[B] 



(y) 




[K] 



(y ' > 

{y ’ ) 



{e ( t) } 




Equations 4.3.6, 4.3.8, and 4.3.11 will be solved in Sections 
6 and 7. 

4.4 Three Degree of Freedom: Rotational. Horizontal, 

and Vertical Translation Response 

4.4.1 Approach 

This section develops a three degree of freedom 
mathematical model of a drydocked submarine. The idealized 
model of the vessel-drydock system is shown in Figure 4.8. In 
this system, vessel motion is restricted to rotation about the 
keel, and in horizontal transverse and vertical translations 
relative to the keel and side pier supports. The approach for 
determining vessel-drydock failure is identical to the two 
degree of freedom case. 

The three degree of freedom mathematical model predicts 
and analyzes submarine rotational, horizontal and vertical 
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transitional responses caused by seismic ground accelerations 
normal to the vessel's longitudinal axis. The rotational 
response is in radians, the translational response in inches, 
and all three responses are expressed as a function of time. 
The keel and side pier forces generated by the responses are 
compared to applicable pier failure criteria found in section 
3 to check blocking stability. 

4.4.2 Three Degree of Freedom Equation of Motion Model 
As in the two degree of freedom case, the keel is the 
origin for the three degree of freedom system with absolute 
rotation, 6, and translation coordinate, x, is the horizontal 
vessel motion relative to the drydock. The linear relative 
translation coordinate, y, is defined with respect to vessel 
and drydock translation relative to the ground. This 
coordinate y describes the vertical motion of the vessel 
hull relative to the drydock and the drydock blocks. Defining 
v as the position of the submarine keel relative to the 
ground, then the following relationships hold; 

y = v - y g 
v = y + y g 

and v=y+y g (4.4.1) 

Likewise, u=x+x g (4.4.2) 

The simplified vessel-block system showing coordinates (x, y 
and 0) at rest and when excited are shown in Figure 4.8 and 
4.9 respectively. 

As developed in Section 4.3.2, the first and third 
equations of motion in the three degree of freedom model are 
MX + MKG0 + c x x + (2K sh + K^x = -Mx g (4.4.3) 
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and 



I k 9 + MKG x - MK5y 9 + C 9 

+ [ (B 2 /2 ) K sv - WKG]0 = -MKG X g + MKG y g 0 . (4.4.4) 

A balance of forces in the y direction, as shown in 
Figure 4.9, yields the second equation of motion 
^ Fy = Mv = (F^sv - Fi dv ) + (F2 S v ” F 2dv) 

+ ( F 3sv " F 3dv> - W. (4.4.5) 



The above equations can be simplified by the following: 

1) The vessel weight, W, must be equal to the summation 
of static block forces so that static equilibrium holds true, 
i.e. 



Fisv + F 2sv F 3SV — W (4.4.6) 

2) The dynamic block forces, F^, as seen from Figure 
4.9 equals to a force due to the modelled springs plus a 
dissipative force due to system damping. These dynamic block 
forces are only a function of vertical displacement, y, and 
velocity, y. The formulation of these forces is identical to 
the procedure used in Section 4.3.2 and will not be repeated. 
The final form of the dynamic block forces, F^, is 

F dv = (2K SV + K kv )y + C y y (4.4.7) 

where K sv = side pier vertical stiffness 
Kkv = keel pier vertical stiffness 
Cy = system vertical damping coefficient. 

Substituting Equations 4.4.1, 4.4.2, 4.4.6 and 4.4.7 into 
Equation 4.4.5 and rearranging components yields 

My + C y y + (2K SV + K kv )y = -My g . (4.4.8) 

The above equation along with Equations 4.4.3 and 4.4.4 



- 40 - 



formulate the three equations of motion for the three degree 
of freedom model . 

The non-linear terms, MKGyQ and MKGyg^, found in 
Equation 4.4.4 will be removed at this time leaving a 
linearized equation since 0 is assumed to be small. The 
effects of the non-linear terms will be evaluated in Section 6 
and 7 . 

Now the system of equations for the three degree of 
freedom is as follows; 

mu x + m 13 0 + x + k^ x = -mu Xg (4.4.9a) 

m 22 y + c 2 y + k 22 y = 22 yg (4.4.9b) 

m 33 0 + m 31 x + c 3 0 + k 33 0 = -m 31 Xg (4.4.9c) 

where = m 22 = M 

m 1 3 — m 3 ^ — MKG 
m 3 3 = I k 

c± = system horizontal damping coefficient 
c 2 = system vertical damping coefficient 
c 3 = system rotational damping coefficient 

k ll = 2K sh + K kh 

k 22 = 2K sv + K kv 
k 33 = B 2 /2 K sv - WKG. 

The first and third equations of motions, Equations 4.4.11a 
and 4.4.11c, respectively, are coupled in the mass times 
acceleration term, which is known as inertial coupling. Also, 
these equations are coupled in the damping term because two of 
the system's three natural modes have horizontal translation 
and rotation involved in them. The second equation of motion, 
Equation 4.4.11b, is uncoupled from the rest of equations. 



- 41 - 



Hence, the last natural mode of the system only depends on 
vertical translation. 

In order to introduce damping, this model uses 5 percent 
of the critical damping coefficient as system damping. To 
evaluate damping, a modal analysis for the system in question 
must be performed. This analysis is performed in Appendix 3 . 
System damping coefficients can be found using this modal 
anlaysis technique, but for now, the coefficients will be left 
in terms of c^, 0^3, C 22 / c 3 i and C33. Substituting these 
coefficients into Equations 4.4.11a, 4.4.11b, and 4.4.11c, the 
equations of motion can be expressed in matrix notation as 
Equation 4.4.12: 

[M] {£_!.} + [C] {y_!_} + [K] {£!} = { E_|_ ( t ) } (4.4.10) 



where 

{y ’ } 



= Response Vector = 



■III 







m ll 


O 


m 13 


[M] 


= Mass Matrix = 


O 


m 22 


0 






m 31 


O 


™33 








k ll 


O 



[K] 



= Stiffness Matrix = 



k 22 0 



{E' (t) } 


= Seismic 


Forcing Vector 


= (" m ll 










j m 22 










(“ m 31 


[C] 


= Damping 


Matrix = 


C 11 


o c 13 








0 


C 22 0 








C 31 


° C 33 



'33_ 

X 

9 

9 9 

x 
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To simplify the explanation the matrices will be redefined as 
[A] <y} + [B] {y> = {E (t) } (4.4.11) 

where 



[A] = 



EB] = 



{y} = 



{ E (t ) } = 



0 


EM] 


[M] 


[C] 


- EM] 


0 


0 


[K] 



\y'\] 

0 

{E' (t) } 



Equations 4.4.3, 4.4.4, 4.4.8, and 4.4.11 will be solved in 
Section 6 and 7. 
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Vessel-Drydock System with Applied Quasi- static Force 



Fs -F h 







FIGURE 4.2 

Vessel-Drydock System Force Diagram with Applied Quasi-static 

Force 



-4*i- 




FIGURE 4.3 

Idealized One Degree of Freedom Vessel- Drydock System at Rest 




Idealized One Degree of Freedom Vessel- Drydock System, 

Excited 
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FIGURE 4.5 

Dynamic Side Block Force 
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FIGURE 4.6 

Idealized Two Degree of Freedom Vessel-Drydock System at Rest 




Xg 

FIGURE 4.7 

Idealized Two Degree of Freedom Vessel-Drydock System, 

Excited 
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FIGURE 4.8 

Idealized Three Degree of Freedom Vessel-Drydock System at 

Rest 




Idealized Three Degree of Freedom Vessel-Drydock System, 

Excited 
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5.0 SYSTEM PARAMETERS 



In order to study the response of a vessel-drydock system 
to seismic loading, eleven typical configurations have been 
selected for implementation into the one, two and three 
degrees of freedom models developed in the previous section. 
This section explores these eleven systems along with the 
modeling of the pier blocks for stiffness. 

5.1 Vessel-Drvdock System Parameters 

Eleven drydock system configurations have been chosen for 
analysis and are defined in Table 5.1. The vessel-drydock 
parameters, corresponding to each system, were taken from the 
appropriate NAVSEA drawings . 

TABLE 5.1 

Vessel-Drydock Configurations 



System 


Hull 


Block 

Type 


Longitudinal 
Block Sprinq 


NAVSEA 

Drawinq 


1 


616 


Composite 


8 ft 


845-2006640 


2 


616 


Composite 


16 ft 


845-2006640 


3 


616 


Timber 


8 ft 


845-2006640 


4 


616 


Timber 


16 ft 


845-2006640 


5 


616 


Timber Side/ 
Concrete Keel 


16 ft 


845-2006640 


6 


726 


Composite 


8 ft 


845-4862749 


7 


726 


Composite 


12 ft 


845-4862749 


8 


726 


Composite 


16 ft 


845-4862749 


9 


688 


Composite/Timber 


12 ft 


845-4403511 


10 


637 


Composite/Timber 


12 ft 


845-2140554 


11 


637 


Compos ite/Timber 


12 ft 


845-2140554 
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Table 5.2 lists the vessel radii, weight, distance from center 
of gravity to keel, KG, moments of inertia about vessel center 
of gravity, I c .g. and keel, 1^, and side block transverse 
spacing, B. 



TABLE 5.2 

Key Vessel Parameters 

Hull 



SYSTEM 


Radius 

(inches) 


Weight 

(Kips) 


KG 

(inches) 


leg ~ 
(K-in-sec 2 ) 


?k ? 

(K-in-sec z ) 


B 

( inches 


1-5 


198 


16,396 


193 


831,257 


2 ,411,000 


144 


6-8 


252 


37,656 


223 


3,097,535 


7,949,000 


180 


9 


198 


13 , 624 


193 


691,852 


2,007,000 


138 


10-11 


190 


9,529 


174 


445,592 


1,193,000 


138 



5.2 Block Parameters 11 



The models developed in Section 4 require the evaluation 
of the block's vertical and horizontal spring stiffness 
constants, K v and K^. The spring stiffness is the amount of 
force required to produce one unit of displacement in the 
respective direction. In order to facilitate the use of 
Hooke's law in evaluating spring constants, the assumption of 
homogeneous isotropic behavior will be used for the block 
materials: wood and concrete. Both the side and keel piers 

must be analyzed for K v and %. 

5.2.1 Vertical Stiffness Ky 11 
Using Hooke's Law, the linear relationship between 
vertical applied force and deformation can be written as 

F _ Eh 
A H 
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where F = vertical force 

A = area under applied force 
E = modulus of elasticity 
h = change in height 
H = original height. 



Hence, the spring constant for a given material is 

, F EA 

v h H ' 



(5.1) 



For the eleven drydock configurations, two types of 
standard blocks (composite and timber) are used. A standard 
composite block which is composed of a softwood cap, a 
hardwood middle portion, and a concrete bottom can be modelled 
as a series of three springs (Figure 5.1). The resultant 
block spring stiffness is 

^v block = t (1/^cap) + (l/^oak) ( V^con) ] ^ (5*2) 

where 



k ca p = softwood cap spring constant 
k oa k = hardwood middle spring constant 
k con = concrete bottom spring constant. 

All three constants listed above are calculated from 
Equation 5.1. A standard timber block is composed of a softwood 
cap and a hardwood body (Figure 5.2). Its block spring 
constant is 

^v block = [ ( Vh C ap) (i/^oak) 3 • (5.3) 

The total vertical stiffness for a drydock configuration 
can be computed by multiplying k v block times the number of 
blocks. The keel pier has a total vetical stiffness of 
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K kv - block * n k 

and for the side pier 

K sv = block * n s 

whre and n s are the number of blocks in one row for their 
respective pier. 

The values of total vertical stiffness for both keel and 
side blocks of the eleven drydock configurations are found in 
Table 5.4. The following parameters are used in the 
formulation of Table 5.4: 

A = cross-sectional area of softwood cap 

E cap = 22.5 ksi 

E oak = 31.675 ksi 

E con = 2,000 ksi 

H cap “ 4 inches 

H oa k = 29 inches for composite block 
= 56 inches for timber block 
H CO n = 27 inches. 

5.2.2 Horizontal Stiffness K^ 11 
In order to determine the horizontal stiffness, two types 
of deformation must be looked at. They are the block-cap 
displacements due to bending and shear deformations. Modelling 
as a continuous cantilever beam subjected to a concentrated 
lateral force at the cap surface (Figure 5.3), the bending 
displacement due to the applied force P of a composite block is 

P(H 1 +H 2 +H 3 ) P(E 1 -E 2 ) (H 2 +H 3 ) 

d B = 3E^I + 3E 1 E 2 I 
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P(E 2 I-E 3 I 3 )H3 

3E 2 E 3 II 3 



(5.4) 



where Hi = height of concrete (27 inches) 

H 2 = height of hardwood (29 inches) 

H 3 = height of softwood cap (4 inches) 

Ei = modulus of elasticity of concrete 

E 2 = modulus of elasticity of hardwood 

E 3 = modulus of elasticity of softwood cap 

I = moment of inertia of block's cross-section 

I 3 = moment of inertia of cap's cross-section. 

For the bending displacement of a timber block, Equation 5.4 
still holds true with the alterations: 

= height of hardwood (27 inches) 

E^ = modulus of elasticity of hardwood. 

In shear, deformation can be determined by modelling the 
composite block as an element subjected to shear stress at the 

top (Figure 5.4). The shear displacement becomes 

(1+v ) PH, 2 (1+v , ) PH- 

, _ con 1 wood 2 

s AE AE_ 



+ 



2(1+v „ood )PH 3 

AE_ 



(5.5) 

where A = area of cap's cross section 

v wood = Poisson's ratio for wood (0.30 is used) 
v con = Poisson's ratio for concrete (0.15 is used). 

By redefining Hi and Ei, Equation 5.5 yields the shear 
displacement of a timber block. 
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The horizontal spring constant of any one block can be 
determined through the Hooke ' s Law relation 

k h = P/(d B + d s ) (5.6) 

The total horizontal stiffnesses for the keel and side piers 
are 

K kh = k h * n k 
and K sh = k h * n s . 

The values for the total area of the rows of blocks for 
each of the eleven systems are listed in Table 5.3. Also listed 
are the corresponding total area of the cap blocks for a row of 
side or keel blocks. The values of total horizontal and 
vertical stiffnesses of the keel and side piers for the eleven 
drydock configurations are listed in Table 5.4. 
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Keel Pier Cross-Sectional Area Analysis 
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TABLE 5.3 (continued) 
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TABLE 5.4 



Keel and Side Pier Stiffness in Pounds Per Inch 



Kkv 


K S v 


Kkh 


K sh 


50089159 


7545899. 


14018741 
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SOFTWOOD CAP 




FIGURE 5.1 

Standard Composite Block Stiffness Model 




FIGURE 5.2 

Standard Timber Block Stiffness Model 
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ir\ 




FIGURE 5.3 

Standard Block in Bending 



i r\ 




FIGURE 5.4 

Standard Block in Shear 
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6.0 VESSEL RESPONSE TO SEISMIC EXCITATION 



This portion of the thesis discusses methods of 
determining the seismic response of the vessel-drydock system 
mathematical models developed in Section 4. The numerical 
analysis schemes to be implemented are a modal analysis method 
(Section 6.1) which looks at the linearized equations of 
motion, and a fourth order Runge-Kutta numerical method 
(Section 6.2) which analyzes the non-linear aspect and also 
verifies the linear solution. Section 6.3 incorporates a 
response spectrum analysis to solve the mathematical models. 
Finally, Section 6.4 generates system response for the eleven 
vessel-drydock configurations using the quasi-static force 
method. 

6 . 1 Vessel Response Using Modal Analysis 

In order to use the modal analysis method described 
herein, only the linearized equations of motion mathematical 
models can be evaluated. These are Equations 4.2.8, 4.3.11 
and 4.4.11 which correspond to the one, two, and three degree 
of freedom models respectively. Since all three linearized 
equations are in the form 

[A]{y} + [B] {y} = {E(t)>, 

the three degree of freedom model will be analyzed. The other 
two cases are just a reduced version of the third. 

To solve the linearized three degree of freedom 
equations of motion (4.4.11) in matrix notation and obtain 
system responses as a function of time, the problem is 
solved as an eigenvalue problem with a procedure described in 
reference [13]. As a first step, the equations are decoupled. 
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To accomplish this, first consider the free vibration problem 



[A]{£} + [B]{y} = 0 

Assuming Equation 6.1.1 has a solution in the form 



( 6 . 1 . 1 ) 



{*} = {4^} e Xt , 



Equation 6.1.1 can be expressed as 
X[A]{^} + [B] { i/j ) = 0. 
where A = complex eigenvalues 
{ 4 1 } = modal vectors. 

Solving the matrix determinant problem 

I MA] + [B] | = 0 



( 6 . 1 . 2 ) 



(6.1.3) 



yields the eigenvalues. In the three degree of freedom case, 
[A] and [B] are 6x6 matrices hence the solution of Equation 
6.1.3 will yield 6 eigenvalues. These eigenvalues will be 
comprised of three sets of complex conjugates. Now the modal 
vectors can be found by substituting the eigenvalues back into 
Equation 6.1.2 and solving. The modal vectors are then 
combined to form the modal matrix [ where 

m = [(Ml (W 2 ( 4 1 ) 3 { W 4 {40 5 {Wei- 
The equations of motion can then be decoupled by using the 
modal matrix , i . e . 



[4o T [B]m = c — b — j 

where [ — -A ] and [ — > B ] are diagonal matrices. 

Once the modal matrix is found, the forced vibration 
problem is analyzed. Assume the following 

iy) = m {z} 

and {z} = [4J]" 1 {y}, 



U] T [A] (40 = [ -A — ] 



and 
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where z is a complex coordinate system. Equation 4.4.1 is now 
expressed as 



[A] [i>] {z} + [B] [i|>] {z} = {E (t ) } . 

Using [i|;] T , the above equation can now be in the form 
r-A^] {z} + [ B — .] {z} = m T {E} — {N (t) } . 

This will yield six uncoupled equation of motion in the form 
a ii Zj_(t) + bj_j_ Zj_(t) = Nj_ (t) , i = 1,6 
or Zi(t) - Zi(t) = (1/aii) Ni(t) (6.1.4) 

where v j_ = 

The solution of Equation 6.1.4 is 



Zi(t) = 



l ii 



/ e vi(t- 



t) 



Once Equation 6.1.5 is solved, 
relationship 



Ni(x) dx, i = 1,6. (6.1.5) 

the {y} matrix is found by the 



{y} = [^] {z}. 

Recalling from Section 4.4.2 the terms of the three 
degree of freedom equations of motion {y} matrix is 

x 



y 

0 



{y> 






y 

0 



Hence, the desired vessel-drydock system response can be found 
in the fourth, fifth and sixth terms of the {y} matrix, i.e. 
x, y and 6. 
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6.2 Vessel Response: Fourth Order Runge-Kutta Numerical 

Method 

A fourth order Runge-Kutta numerical scheme is chosen to 
analyze the non-linear equations of motion due to its 
simplicity and easy use of computer programming. It also will 
be used to verify solutions found by other methods. The 
non-linear Equations 4.2.6; 4.3.6 and 4.3.8; 4.4.3. 4.4.4 and 
4.4.8 correspond to the one, two and three degree of 
freedom models respectively. Since all three systems of 
equations are similar in form, only the three degree of 
freedom non-linear model will be developed. 

The three degree of freedom non-linear equations of 
motion are 



• • • 
m^i x + mj_ 3 6 + c^i x + Cj _3 6 + x — — m^j Xg 

m 22 y + c 22 y + k 22 Y = ~ m22 V< 3 



and 



*33 6 + m 31 x - m 31 y 9 + c 33 6 + c 31 x + k 31 e 

= -m 31 x g + m 31 y g e . 

The coefficient of the above equations can be found in Section 



4.4.2. These equations can be rearranged into 



m 13 .. _ _ c ll • _ ^13 • _ ^11 _ 

X + m, , 9 m 11 X m 11 6 g 



11 



c y 

c 22 • 22 

y v y - y„ 

^ m 22 ‘ ra 22 9 



(6.1.2a) 



(6.1.2b) 



and 



m. 



. 31 " 

6 + x 



m 



33 



C 3 3 * 

— ±± e 

m 33 
m 31 - 
m 33 9 



C 31 * 
x 

m 33 



-( 



'33 



m 



33 



+ m 31 C 22 
m 33 m 22 



+ 



m 31 k 22 

m 22 m 33 



y)e 



(6.2.1c) 
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